Computational Models for Clinical Applications in Personalized Medicine—Guidelines and Recommendations for Data Integration and Model Validation

The future development of personalized medicine depends on a vast exchange of data from different sources, as well as harmonized integrative analysis of large-scale clinical health and sample data. Computational-modelling approaches play a key role in the analysis of the underlying molecular processes and pathways that characterize human biology, but they also lead to a more profound understanding of the mechanisms and factors that drive diseases; hence, they allow personalized treatment strategies that are guided by central clinical questions. However, despite the growing popularity of computational-modelling approaches in different stakeholder communities, there are still many hurdles to overcome for their clinical routine implementation in the future. Especially the integration of heterogeneous data from multiple sources and types are challenging tasks that require clear guidelines that also have to comply with high ethical and legal standards. Here, we discuss the most relevant computational models for personalized medicine in detail that can be considered as best-practice guidelines for application in clinical care. We define specific challenges and provide applicable guidelines and recommendations for study design, data acquisition, and operation as well as for model validation and clinical translation and other research areas

Investigating the genetic association between ERAP1 and ankylosing spondylitis

A strong association between ERAP1 and ankylosing spondylitis (AS) was recently identified by the Wellcome Trust Case Control Consortium and the Australo-Anglo-American Spondylitis Consortium (WTCCC-TASC) study. ERAP1 is highly polymorphic with strong linkage disequilibrium evident across the gene. We therefore conducted a series of experiments to try to identify the primary genetic association(s) with ERAP1. We replicated the original associations in an independent set of 730 patients and 1021 controls, resequenced ERAP1 to define the full extent of coding polymorphisms and tested all variants in additional association studies. The genetic association with ERAP1 was independently confirmed; the strongest association was with rs30187 in the replication set (P = 3.4 × 10−3). When the data were combined with the original WTCCC-TASC study the strongest association was with rs27044 (P = 1.1 × 10−9). We identified 33 sequence polymorphisms in ERAP1, including three novel and eight known non-synonymous polymorphisms. We report several new associations between AS and polymorphisms distributed across ERAP1 from the extended case–control study, the most significant of which was with rs27434 (P = 4.7 × 10−7). Regression analysis failed to identify a primary association clearly; we therefore used data from HapMap to impute genotypes for an additional 205 non-coding SNPs located within and adjacent to ERAP1. A number of highly significant associations (P < 5 × 10−9) were identified in regulatory sequences which are good candidates for causing susceptibility to AS, possibly by regulating ERAP1 expression

Large-Scale Genome-Wide Meta-Analysis of Polycystic Ovary Syndrome Suggests Shared Genetic Architecture for Different Diagnosis Criteria

Polycystic ovary syndrome (PCOS) is a disorder characterized by hyperandrogenism, ovulatory dysfunction and polycystic ovarian morphology. Affected women frequently have metabolic disturbances including insulin resistance and dysregulation of glucose homeostasis. PCOS is diagnosed with two different sets of diagnostic criteria, resulting in a phenotypic spectrum of PCOS cases. The genetic similarities between cases diagnosed based on the two criteria have been largely unknown. Previous studies in Chinese and European subjects have identified 16 loci associated with risk of PCOS. We report a fixed-effect, inverse-weighted-variance meta-analysis from 10,074 PCOS cases and 103,164 controls of European ancestry and characterisation of PCOS related traits. We identified 3 novel loci (near PLGRKT, ZBTB16 and MAPRE1), and provide replication of 11 previously reported loci. Only one locus differed significantly in its association by diagnostic criteria; otherwise the genetic architecture was similar between PCOS diagnosed by self-report and PCOS diagnosed by NIH or non-NIH Rotterdam criteria across common variants at 13 loci. Identified variants were associated with hyperandrogenism, gonadotropin regulation and testosterone levels in affected women. Linkage disequilibrium score regression analysis revealed genetic correlations with obesity, fasting insulin, type 2 diabetes, lipid levels and coronary artery disease, indicating shared genetic architecture between metabolic traits and PCOS. Mendelian randomization analyses suggested variants associated with body mass index, fasting insulin, menopause timing, depression and male-pattern balding play a causal role in PCOS. The data thus demonstrate 3 novel loci associated with PCOS and similar genetic architecture for all diagnostic criteria. The data also provide the first genetic evidence for a male phenotype for PCOS and a causal link to depression, a previously hypothesized comorbid disease. Thus, the genetics provide a comprehensive view of PCOS that encompasses multiple diagnostic criteria, gender, reproductive potential and mental health

Protein-coding variants implicate novel genes related to lipid homeostasis contributing to body-fat distribution.

Body-fat distribution is a risk factor for adverse cardiovascular health consequences. We analyzed the association of body-fat distribution, assessed by waist-to-hip ratio adjusted for body mass index, with 228,985 predicted coding and splice site variants available on exome arrays in up to 344,369 individuals from five major ancestries (discovery) and 132,177 European-ancestry individuals (validation). We identified 15 common (minor allele frequency, MAF ≥5%) and nine low-frequency or rare (MAF &lt;5%) coding novel variants. Pathway/gene set enrichment analyses identified lipid particle, adiponectin, abnormal white adipose tissue physiology and bone development and morphology as important contributors to fat distribution, while cross-trait associations highlight cardiometabolic traits. In functional follow-up analyses, specifically in Drosophila RNAi-knockdowns, we observed a significant increase in the total body triglyceride levels for two genes (DNAH10 and PLXND1). We implicate novel genes in fat distribution, stressing the importance of interrogating low-frequency and protein-coding variants

A polygenic and phenotypic risk prediction for polycystic ovary syndrome evaluated by phenomewide association studies

Context: As many as 75% of patients with polycystic ovary syndrome (PCOS) are estimated tobe unidentified in clinical practice. Objective: Utilizing polygenic risk prediction, we aim to identify the phenome-widecomorbidity patterns characteristic of PCOS to improve accurate diagnosis and preventivetreatment.Design, Patients, and Methods: Leveraging the electronic health records (EHRs) of 124 852individuals, we developed a PCOS risk prediction algorithm by combining polygenic risk scores(PRS) with PCOS component phenotypes into a polygenic and phenotypic risk score (PPRS). Weevaluated its predictive capability across different ancestries and perform a PRS-based phenomewide association study (PheWAS) to assess the phenomic expression of the heightened risk ofPCOS.Results: The integrated polygenic prediction improved the average performance (pseudo-R2)for PCOS detection by 0.228 (61.5-fold), 0.224 (58.8-fold), 0.211 (57.0-fold) over the null modelacross European, African, and multi-ancestry participants respectively. The subsequent PRSpowered PheWAS identified a high level of shared biology between PCOS and a range ofmetabolic and endocrine outcomes, especially with obesity and diabetes: "morbid obesity","type 2 diabetes", "hypercholesterolemia", "disorders of lipid metabolism", "hypertension",and "sleep apnea" reaching phenome-wide significance.Conclusions: Our study has expanded the methodological utility of PRS in patient stratificationand risk prediction, especially in a multifactorial condition like PCOS, across different geneticorigins. By utilizing the individual genome-phenome data available from the EHR, our approachalso demonstrates that polygenic prediction by PRS can provide valuable opportunities todiscover the pleiotropic phenomic network associated with PCOS pathogenesis.Abbreviations: AA, African ancestry; ANOVA, analysis of variance; BMI, body mass index; EA,European ancestry; EHR, electronic health records; eMERGE, electronic Medical Records andGenomics Network; GWAS, genome-wide association study; IBD, identity-by-descent; ICDCM, International Classification of Diseases, Clinical Modification; LD, linkage disequilibrium;MA, multi-ancestry; MAF, minor allele frequency; NIH, National Institutes of Health; PCA,principal component analysis; PheWAS, phenome-wide association study; PCOS, polycysticovary syndrome; PPRS, polygenic and phenotypic risk score; PRS, polygenic risk sc

Exome-Derived Adiponectin-Associated Variants Implicate Obesity and Lipid Biology

Circulating levels of adiponectin, an adipocyte-secreted protein associated with cardiovascular and metabolic risk, are highly heritable. To gain insights into the biology that regulates adiponectin levels, we performed an exome array meta-analysis of 265,780 genetic variants in 67,739 individuals of European, Hispanic, African American, and East Asian ancestry. We identified 20 loci associated with adiponectin, including 11 that had been reported previously (p .60) spanning as much as 900 kb. To identify potential genes and mechanisms through which the previously unreported association signals act to affect adiponectin levels, we assessed cross-trait associations, expression quantitative trait loci in subcutaneous adipose, and biological pathways of nearby genes. Eight of the nine loci were also associated (p <1 x 10(-4)) with at least one obesity or lipid trait. Candidate genes include PRKAR2A, PTH1R, and HDAC9, which have been suggested to play roles in adipocyte differentiation or bone marrow adipose tissue. Taken together, these findings provide further insights into the processes that influence circulating adiponectin levels.Peer reviewe

Genetic Studies of Leptin Concentrations Implicate Leptin in the Regulation of Early Adiposity.

Leptin influences food intake by informing the brain about the status of body fat stores. Rare LEP mutations associated with congenital leptin deficiency cause severe early-onset obesity that can be mitigated by administering leptin. However, the role of genetic regulation of leptin in polygenic obesity remains poorly understood. We performed an exome-based analysis in up to 57,232 individuals of diverse ancestries to identify genetic variants that influence adiposity-adjusted leptin concentrations. We identify five novel variants, including four missense variants, in LEP, ZNF800, KLHL31, and ACTL9, and one intergenic variant near KLF14. The missense variant Val94Met (rs17151919) in LEP was common in individuals of African ancestry only, and its association with lower leptin concentrations was specific to this ancestry (P = 2 × 10-16, n = 3,901). Using in vitro analyses, we show that the Met94 allele decreases leptin secretion. We also show that the Met94 allele is associated with higher BMI in young African-ancestry children but not in adults, suggesting that leptin regulates early adiposity

Oral Abstracts 7: RA ClinicalO37. Long-Term Outcomes of Early RA Patients Initiated with Adalimumab Plus Methotrexate Compared with Methotrexate Alone Following a Targeted Treatment Approach

Background: This analysis assessed, on a group level, whether there is a long-term advantage for early RA patients treated with adalimumab (ADA) + MTX vs those initially treated with placebo (PBO) + MTX who either responded to therapy or added ADA following inadequate response (IR). Methods: OPTIMA was a 78- week, randomized, controlled trial of ADA + MTX vs PBO + MTX in MTX-naïve early (<1 year) RA patients. Therapy was adjusted at week 26: ADA + MTX-responders (R) who achieved DAS28 (CRP) <3.2 at weeks 22 and 26 (Period 1, P1) were re-randomized to withdraw or continue ADA and PBO + MTX-R continued randomized therapy for 52 weeks (P2); IR-patients received open-label (OL) ADA + MTX during P2. This post hoc analysis evaluated the proportion of patients at week 78 with DAS28 (CRP) <3.2, HAQ-DI <0.5, and/or ΔmTSS ≤0.5 by initial treatment. To account for patients who withdrew ADA during P2, an equivalent proportion of R was imputed from ADA + MTX-R patients. Results: At week 26, significantly more patients had low disease activity, normal function, and/or no radiographic progression with ADA + MTX vs PBO + MTX (Table 1). Differences in clinical and functional outcomes disappeared following additional treatment, when PBO + MTX-IR (n = 348/460) switched to OL ADA + MTX. Addition of OL ADA slowed radiographic progression, but more patients who received ADA + MTX from baseline had no radiographic progression at week 78 than patients who received initial PBO + MTX. Conclusions: Early RA patients treated with PBO + MTX achieved comparable long-term clinical and functional outcomes on a group level as those who began ADA + MTX, but only when therapy was optimized by the addition of ADA in PBO + MTX-IR. Still, ADA + MTX therapy conferred a radiographic benefit although the difference did not appear to translate to an additional functional benefit. Disclosures: P.E., AbbVie, Merck, Pfizer, UCB, Roche, BMS—Provided Expert Advice, Undertaken Trials, AbbVie—AbbVie sponsored the study, contributed to its design, and participated in the collection, analysis, and interpretation of the data, and in the writing, reviewing, and approval of the final version. R.F., AbbVie, Pfizer, Merck, Roche, UCB, Celgene, Amgen, AstraZeneca, BMS, Janssen, Lilly, Novartis—Research Grants, Consultation Fees. S.F., AbbVie—Employee, Stocks. A.K., AbbVie, Amgen, AstraZeneca, BMS, Celgene, Centocor-Janssen, Pfizer, Roche, UCB—Research Grants, Consultation Fees. H.K., AbbVie—Employee, Stocks. S.R., AbbVie—Employee, Stocks. J.S., AbbVie, Amgen, AstraZeneca, BMS, Celgene, Centocor-Janssen, GlaxoSmithKline, Lilly, Pfizer (Wyeth), MSD (Schering-Plough), Novo-Nordisk, Roche, Sandoz, UCB—Research Grants, Consultation Fees. R.V., AbbVie, BMS, GlaxoSmithKline, Human Genome Sciences, Merck, Pfizer, Roche, UCB Pharma—Consultation Fees, Research Support. Table 1.Week 78 clinical, functional, and radiographic outcomes in patients who received continued ADA + MTX vs those who continued PBO + MTX or added open-label ADA following an inadequate response ADA + MTX, n/N (%)a PBO + MTX, n/N (%)b Outcome Week 26 Week 52 Week 78 Week 26 Week 52 Week 78 DAS28 (CRP) <3.2 246/466 (53) 304/465 (65) 303/465 (65) 139/460 (30)*** 284/460 (62) 300/460 (65) HAQ-DI <0.5 211/466 (45) 220/466 (47) 224/466 (48) 150/460 (33)*** 203/460 (44) 208/460 (45) ΔmTSS ≤0.5 402/462 (87) 379/445 (86) 382/443 (86) 330/459 (72)*** 318/440 (72)*** 318/440 (72)*** DAS28 (CRP) <3.2 + ΔmTSS ≤0.5 216/462 (47) 260/443 (59) 266/443 (60) 112/459 (24)*** 196/440 (45) 211/440 (48)*** DAS28 (CRP) <3.2 + HAQ-DI <0.5 + ΔmTSS ≤0.5 146/462 (32) 168/443 (38) 174/443 (39) 82/459 (18)*** 120/440 (27)*** 135/440 (31)** aIncludes patients from the ADA Continuation (n = 105) and OL ADA Carry On (n = 259) arms, as well as the proportional equivalent number of responders from the ADA Withdrawal arm (n = 102). bIncludes patients from the MTX Continuation (n = 112) and Rescue ADA (n = 348) arms. Last observation carried forward: DAS28 (CRP) and HAQ-DI; Multiple imputations: ΔmTSS. ***P < 0.001 and **iP < 0.01, respectively, for differences between initial treatments from chi-squar